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ABSTRACT 

New photometric solutions have been carried out on the important eccentric eclipsing system 
V380 Cygni (B1.5 II-III + B2 V) from UBV differential photoelectric photometry obtained by 
us. The photometric elements obtained from the analysis of the light curves have been combined 
with the spectroscopic solution recently published by Popper & Guinan and have led to the 
physical properties of the system components. The effective temperature of the stars has been 
determined by fitting IUE UV spectrophotometry to Kurucz model atmospheres and compared 
with other determinations from broad-band and intermediate-band standard photometry The 
values of mass, absolute radius, and effective temperature, for the primary and secondary stars 
are: 11.1±0.5 M , 14.7±0.2 R , 21 350±400K, and 6.95±0.25 M Q , 3.74±0.07 R Q , 20 500±500 K, 
respectively. In addition, a re-determination of the system's apsidal motion rate has been done 
from the analysis of 12 eclipse timings obtained from 1923 to 1995. The apsidal motion study 
yields the internal mass distribution of the more luminous component. Using stellar structure 
and evolutionary models with modern input physics, tests on the extent of convection in the core 
of the more massive B1.5 II-III star of the system have been carried out. Both the analysis of 
the log g — log T c g diagram and the apsidal motion study indicate a star with a larger convective 
core, and thus more centrally condensed, than currently assumed. This has been quantified in 
form of an overshooting parameter with a value of a ov ~ 0.6 ± 0.1. Finally, the tidal evolution 
of the system (synchronization and circularization times) has also been studied. 

Subject headings: stars: fundamental parameters — binaries: eclipsing — stars: evolution — stars: 
early-type — stars: individual (V380 Cyg) 
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1. Introduction 

The bright eclipsing binary V380 Cyg (HR 7567; 
HD 187879; HIP 97634; V max = 5.68; P = 
12.426 days; B1.5 II-III + B2 V; Hill & Batten 
1984) has properties that make it an important 
"astrophysical laboratory" for studying the struc- 
ture and evolution of massive stars. In particular, 
the extent of convection in the stellar core and the 
internal mass distribution of the primary com- 
ponent can be probed because the system has 
an eccentric orbit with a well-established apsi- 
dal motion rate. Also, V380 Cyg can provide 
independent measures of the initial fractional he- 
lium abundance of the system (Y), which is an 
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important and fundamental quantity but empiri- 
cally difficult to measure. Accurate fundamental 
physical properties of the components (e.g., mass, 
radius, effective temperature, etc.) are however 
required to carry out such analyses. 

V380 Cyg consists of an evolved, more mas- 
sive, and more luminous primary component and 
a main sequence secondary star. It has an ec- 
centric orbit (e = 0.23) and an orbital period of 
12.426 days. Several spectroscopic studies have 
been carried out to determine its orbital proper- 
ties and the masses, temperatures, and luminosi- 
ties of the component stars (see Batten 1962; Pop- 
per 1981; Hill & Batten 1984; Lyubimkov et al. 
1996). The first light curve of the system by Kron 
(1935), as well as the more modern photometry 
by Semeniuk (1968) and Battistini, Bonifazi, & 
Guarnieri (1974), revealed shallow eclipses with 
depths of ~ 0.12 mag and ~ 0.09 mag for pri- 
mary and secondary minima, respectively. More- 
over, these photometric measurements show that 
the secondary eclipse is displaced from 0?5, indi- 
cating an eccentric orbit, while changes in the dis- 
placement between primary and secondary min- 
ima indicate the presence of apsidal motion with 
a period of about 1500 yrs. 

As discussed by Gimenez (1984) and Gimenez, 
Claret, & Guinan (1994), V380 Cyg is an ideal bi- 
nary system for the study of convective overshoot- 
ing in the cores of massive stars because of the evo- 
lutionary stages of its component stars. Because 
of the eccentric orbit and the eclipsing nature of 
the binary, it is possible to determine the apsidal 
motion rate. From this, additional constraints on 
the internal mass distribution and the evolution- 
ary state of an evolved massive star also can be 
established. 

Because V380 Cyg provides potentially impor- 
tant tests of stellar structure and evolution, we 
have carried out new high signal-to-noise spectro- 
scopic and photometric observations of the system 
and performed a new, detailed investigation of its 
properties. In §2 we present the new observations. 
The effective temperature determination, based on 
both standard photometry and UV/optical spec- 
trophotometry, is discussed in §3. In §4 we con- 
centrate on the analysis of the light curves. §5 is 
devoted to the study of eclipse timings, leading to 
an accurate determination of the apsidal motion 
rate. The properties of the system components 



are compared with the predictions of theoretical 
models in §6. The tidal evolution of the system 
(circularization and synchronization times) is an- 
alyzed in §7. Finally, the main conclusions of our 
study are presented in §8. 

2. Observational Data 

New high-quality spectroscopic and photomet- 
ric observations have been collected for V380 Cyg. 
The following sections contain a brief description 
of these data. 

2.1. Spectroscopy 

V380 Cyg is difficult to resolve spectroscop- 
ically into a double-line binary because of the 
large luminosity difference between the compo- 
nent stars. Therefore, very high signal-to-noisc 
and high-resolution spectroscopic observations are 
needed to measure the radial velocity shifts of the 
faint lines of the secondary star. Such high signal- 
to-noise and high-resolution spectroscopy was re- 
cently secured by us and led to a new, accurate 
radial velocity curve, which was carefully analyzed 
in a previous paper (Popper & Guinan 1998; here- 
after Paper I). The spectroscopic solution obtained 
in Paper I (basically velocity semi-amplitudes and 
systemic velocity) has been adopted in the present 
study. 

2.2. Photometry 

UBV differential photoelectric photometry of 
V380 Cyg was carried out from 1988 to 1997 
using robotic, automatic photoelectric telescopes 
(APTs) located on Mt. Hopkins, Arizona. From 
1988 to 1989 the photometry was conducted with 
the Phoenix-10 APT telescope, a 25-cm robotic 
telescope equipped with a 1-P21 photometer and 
filters matched to the standard UBV system. 
During 1995-1997 the photometry was conducted 
using the Four College Consortium (FCC) 0.8- 
m APT. The FCC APT is equipped with filters 
closely matching the standard UBV system and 
a refrigerated Hamamatsu photoelectric detector. 
Differential photometry was obtained in which 
HD 189178 (B5 V, V = 5.47) was employed as the 
comparison star, while HD 188892 (B5 IV, V = 
4.94) served as the check star. No evidence of sig- 
nificant light variations (down to the few milimag- 
nitude level) was found for the comparison-check 
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star sets. The observations were reduced using 
photometric reduction programs at Villanova Uni- 
versity. Differential extinction corrections were 
applied, but these corrections were typically very 
small since the comparison star is within 2 deg. 
of V380 Cyg and the stars were never observed at 
high air mass. The photometric measurements are 
listed in Table 1. 

U, B, and V differential light curves are 
shown in Figure 1. We adopted the follow- 
ing cphcmcris for computing the orbital phase: 
T Mi „i = HJD2441256.544 + 12.425719 E (sec §5). 
As can be seen in the plots, the shallow sec- 
ondary minimum occurs at 0?403. Also, the light 
curves show out-of-eclipse brightness variations 
with maximum light occurring near 0?15. The 
displacement of the secondary eclipse from phase 
0.5, as noted in §1 above, is caused by the eccen- 
tric orbit of the system, while the outside eclipse 
light variations arise primarily from the changing 
tidal distortion, chiefly from the larger star as it 
moves in the eccentric orbit. According to the re- 
cent spectroscopic solution of Paper I, periastron 
passage occurs near 0?15 where the stars reach 
their greatest tidal deformations. 

3. Temperature Determination 

A complete analysis of the V380 Cyg system re- 
quires precise measurements of the stellar effective 
temperatures (T c g). For this purpose, we explored 
the use of published T c ff vs. intrinsic color cali- 
brations utilizing several different photometric 
systems — and also we performed a detailed study 
of the shape of the UV/optical stellar energy dis- 
tribution. Both these methods for determining ef- 
fective temperatures are discussed below. 

Intermediate-band Stromgrcn and Vilnius pho- 
tometric indices for V380 Cyg are available from 
the literature and are listed in Table 2 along with 
our new measurements in the Johnson UBV sys- 
tem. The Johnson photometry was secured out- 
of-eclipse, when both stars were fully visible. The 
phases at which the Stromgren and Vilnius mea- 
surements were made is not known, although this 
is not likely to be a critical problem for V380 Cyg 
because the secondary component is much less lu- 
minous than the primary component and the tem- 
peratures and color indices of the two stars are 
similar. A summary of the reddening and T c g re- 



sults for the three photometric systems is given in 
Table 2. An a priori estimate of the expected un- 
certainties in these results is difficult, but clearly 
depends on both random errors in the stellar pho- 
tometry and systematic errors in the various cali- 
brations themselves. It is clear, however, that the 
level of agreement among the various systems is 
unsatisfactory, with a range of over 4000 K and a 
standard deviation of ^1800 K — nearly 10% of 
the value of T e s itself. 

In view of the discrepancies among the pho- 
tometric results for T ff, we employed a different 
technique for determining this important quantity. 
We modeled the observed shape of the UV/optical 
energy distribution of V380 Cyg using theoreti- 
cal spectra generated with the Kurucz's ATLAS9 
model atmosphere program. This analysis was 
developed by Fitzpatrick & Massa (1999; here- 
after FM99) and has been applied previously to an 
eclipsing binary system by Guinan et al. (1998). 

For a binary system, the observed energy dis- 
tribution depends on the surface fluxes of the bi- 
nary's components and on the attenuating effects 
of distance and interstellar extinction. This rela- 
tionship can be expressed as: 

/Affi = | F P + ^i!sj 2 F sj xl0 -0.4E(B-V)[fc(A-V) + H(V)] 

(1) 

where and are the surface fluxes of the 
primary (P) and secondary (S) components, Rp 
and i?s are the absolute radii, and d is the distance 
to the binary. The last term carries the extinction 
information, including E(B — V), the normalized 
extinction curve k(X - V) = E(X - V)/E(B - V), 
and the ratio of selective-to-total extinction in the 
V band R(V) = A(V)/E(B - V). 

The modeling process consists of a non-linear 
least squares determination of the optimal val- 
ues of all the parameters which contribute to the 
right side of equation 1. In principle, the problem 
can involve solving for two sets of Kurucz model 
parameters (i.e., one set of T c g, surface gravity, 
metallicity [m/H], and microturbulcncc velocity 
fmicro values for each star), the ratios (Rp/d) 2 
and R s /Rp, E(B - V), k(X - V), and R(V). For 
V380 Cyg, however, a number of simplifications 
can be made: (1) the temperature ratio of the 
two stars is known from the light curve analysis 
(see §4 below); (2) the surface gravities are known 
from the combined result of the light and radial 
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velocity curve analyses (see §4); (3) [m/H] can be 
assumed to be identical for both components; (4) 
Wmicro can be assumed to be for the much fainter 
secondary star; (5) the ratio Rg/ Rp is known from 
the light curve analysis (sec §4); (6) the proper- 
ties of normalized UV/optical extinction are con- 
strained by a small number of parameters (Fitz- 
patrick & Massa 1990; Fitzpatrick 1999); and (7) 
the standard value of R(V) = 3.1 may be assumed, 
without affecting the T e g determination. 

Adopting these simplifications, we modeled 
the observed UV/optical energy distribution of 
V380 Cyg by solving for the best-fitting values 
of T P , [m/H] PS , < icro , (R P /d)\ E(B - V), and 
k(X — V). The observational data used in the 
analysis includes the Johnson and Stromgren op- 
tical photometry listed in Table 2 and UV spec- 
trophotometry from the International Ultravio- 
let Explorer (IUE) satellite. The IUE data con- 
sist of a low-dispersion short-wavelength spec- 
trum (SWP38553; 1150 A to 1980 A) and a low- 
dispersion long-wavelength spectrum (LWP17712; 
1980 A to 3000 A), both secured out-of-eclipse. 
These data were obtained from the IUE Fi- 
nal Archive and thus were processed using the 
NEWSIPS system (Nichols & Linksy 1996). We 
utilized the algorithms of Massa & Fitzpatrick 
(2000) to correct these NEWSIPS data for strong 
thermal and temporal systematics and to place 
them on the Hubble Space Telescope absolute 
UV/optical calibration system (Bohlin 1996). The 
final IUE spectrum was then re-samplcd to match 
the wavelength binning of the ATLAS9 models 
(see FM99). This binned spectrum can be seen 
in Figure 2 (small filled circles). The presence of 
the 2175 A interstellar medium extinction feature 
is clear, as are the strong stellar absorption fea- 
tures in the 1400 A, 1500 A, and 1900 A regions 
characteristic of evolved early B-type stars. For 
display purposes, we also show the Johnson and 
Stromgren magnitudes, converted to fluxes, in Fig- 
ure 2 (large filled circles and triangles). However, 
in the fitting procedure the photometric indices 
themselves (i.e., V, B — V, U — B, b — y, mi, c\, 
j3) were compared with synthetic values for the 
models. The calibration of the synthetic photom- 
etry is discussed by FM99 and in a forthcoming 
paper by E. L. Fitzpatrick & D. Massa (in prepa- 
ration). 

The best-fitting, reddened, and distance-at- 



tenuated ATLAS9 model of V380 Cyg's energy 
distribution is shown in Figure 2 (histogram-style 
curve) overplotted on the V380 Cyg observations. 
In the UV region, the residuals to the fit ("Model 
- Star" ) are shown below the energy distribution, 
with the vertical bars indicating the statistical 
error in the IUE data. The major discrepancies 
occur near the interstellar Ly a line (1215 A) and 
the C IV stellar wind line (1550 A). These re- 
gions were excluded from the fit, as indicated by 
the crosses over the data points. Several other 
data points, corresponding to spurious feature in 
the models were also excluded (see FM99). (Note 
that the fit appears poorest in the 2000 to 2300 A 
region because of the low signal-to-noise level of 
the data.) 

The parameters characterizing the fit are listed 
in Table 3. The values of T off and E(B - V) are 
within the range found from the photometric anal- 
ysis (see Table 2) but are much more precisely 
determined. (The synthetic photometry of the 
best-fitting model reproduces all the Johnson and 
Stromgren photometric indices to within the ex- 
pected errors of the data.) The attenuation factor 
of the energy distribution is also derived from the 
analysis, and it yields the distance to the system 
(Table 4) when combined with the absolute radius 
of the primary component (see Guinan et al. 1998 
for further discussion on this issue). The values 
of [m/H] and w m i cr o bear some comment. FM99 
found that in B-type stars with significant stellar 
winds, correlated increases in the best-fitting val- 
ues of Vmicro and decreases in [m/H] occur. These 
combinations reproduce the observed opacity fea- 
tures very well, but it is likely that neither pa- 
rameter faithfully represents a corresponding stel- 
lar property. This probably arises because the 
stellar opacity features form in a region with a 
weak velocity gradient, while the models explic- 
itly assume an atmosphere in hydrostatic equi- 
librium. (Sec FM99 for additional discussion.) 
The implication for this work is that the value 
of [m/H] derived from this analysis is almost cer- 
tainly an underestimate of the true metallicity of 
V380 Cyg. For further analysis, we adopt a metal 
abundance of [m/H] = —0.2 based on recent anal- 
yses of nearby B-type stars (Kilian, Montenbruck, 
& Nissen 1994), solar-type stars (Edvardsson ct 
al. 1993) and interstellar medium studies (Meyer, 
Jura, & Cardelli 1998). 
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Note that this fitting procedure relied on results 
generated by the light and radial velocity curve 
analyses described below, while these other anal- 
yses rely (albeit very weakly) on the T e g's deter- 
mined here. We thus performed these three anal- 
yses in an iterative fashion until a self-consistent 
set of parameters resulted. 

4. Modeling the Light Curves 

The light curves were solved using an improved 
version of the Wilson & Devinney (1971; here- 
after WD) program (Milone, Stagg, & Kurucz 
1992; Milone et al. 1994) that includes Kurucz 
ATLAS9 atmosphere models for the computation 
of the stellar radiative parameters. A detached 
configuration (as suggested by the eccentric or- 
bit) with coupling between luminosity and tem- 
perature was chosen when running the solutions. 
Both reflection and proximity effects (i.e., tidal 
distortion) were taken into account, because the 
light curves clearly show their significance. The 
bolometric albedo and the gravity brightening co- 
efficients were set to a value of 1.0 as usually 
adopted for the radiative envelopes of these early- 
type stars. The mass ratio (q = Ms/Mp) was 
fixed to the spectroscopic value of q = 0.626 and 
the temperature of the primary star was set to 
21 350 K. Finally, pseudo-synchronization (syn- 
chronization at periastron) was assumed for both 
components (see Sect. 7). 

All three UBV light curves were solved simul- 
taneously to achieve a single, mutually consistent, 
solution. In this simultaneous solution, the rel- 
ative weight of each curve, which has significant 
importance, was set according to the r.m.s. resid- 
ual of the comparison— check star sets (i.e., 0.015, 
0.006, and 0.006 for UBV, respectively). The free 
parameters in the light curve fitting were: the ec- 
centricity (e), the longitude of the periastron (u>), 
the phase offset (4> ), the orbital inclination (i), 
the temperature of the secondary (Tg), the grav- 
itational potentials (Op and fig), and the lumi- 
nosity of the primary (Lp). Depending on the 
initial conditions, two different sets of parame- 
ters providing good fits to the observations were 
found. There is not sufficient information in the 
light or radial velocity curves to discriminate be- 
tween these two possible solutions, and additional 
data sources have to be considered to break the 



degeneracy. Since the luminosity ratio between 
the components predicted by the two solutions is 
very different (L§/Lp « 0.05 and L$/Lp « 0.5), 
the spectroscopic light ratio (Lg/Lp « 0.05, from 
Paper I) was used to confirm the first solution as 
the physically preferable one. 

An automatic procedure based on midas scripts 
and FORTRAN programs was employed to execute 
the WD code. At each iteration, the differential 
corrections were applied to the input parameters 
to build the new set of parameters for the next it- 
eration. The result was carefully checked in order 
to avoid possible unphysical situations (e.g., Roche 
lobe filling in detached configuration). A solution 
was defined as the set of parameters for which the 
differential corrections suggested by the WD code 
were smaller than the standard errors during three 
consecutive iterations. However, when a solution 
was found, the program did not stop. Instead, 
it was kept running for a given number of itera- 
tions, 100 in our case. This was done to test the 
stability of the solutions, to evaluate their scat- 
ter, and to check for possible spurious solutions. 
For V380 Cyg, the program rapidly converged to a 
very stable minimum and 47 individual solutions 
resulted from the 100 iterations. A summary of 
the best-fitting parameters to the light curve is 
presented in Table 3. 

We also considered - and ultimately rejected 
- the possibility of a third light contribution to 
the system. When searching both the SIMBAD 
database and the Tycho catalogue (ESA 1997) no 
objects were found within a radius of 1 arcminutc 
centered at the position of V380 Cyg. The system 
was observed by the Hipparcos mission and there 
is no indication of visual duplicity down to the 
spatial and intensity resolution of the satellite (p « 
0.1 — 0.15 mas and Amag « 3.5 — 4). Moreover, no 
evidence of stationary spectral lines belonging to 
a third component are observed in the high S/N 
spectra used for the spectroscopic study. Since the 
contribution of the secondary component (whose 
lines are seen in the spectra) to the total flux of 
the system in the considered wavelengths is only 
about 6 percent, we therefore conclude that any 
third light can be safely neglected. 

Additional WD runs were performed starting 
from the same initial conditions but fixing Tp to 
its proper value ±la. This showed that our re- 
sults from the light curve solution are insensitive 
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to temperature variations within this range. 

A remark should be made at this point concern- 
ing the physical meaning of the primary and the 
secondary eclipses. Since the massive and larger 
component is also somewhat hotter (and so it has 
slightly higher surface flux), one would therefore 
expect the primary (and deeper) eclipse to be a 
transit of the smaller star in front of the bigger 
one. However, this is not the case for V380 Cyg, 
and the primary eclipse turns out to be an occul- 
tation of the smaller and cooler component. Such 
a behavior can be explained in terms of the ec- 
centric orbit and the variations induced by the 
proximity effects, which can be seen clearly as an 
out-of-eclipse wave-like variation in Figure 1. In- 
deed, eccentricity-induced proximity effects pro- 
duce minimum flux near the phase of the transit, 
making it the primary eclipse. 

The uncertainties in the parameters presented 
in Table 3 were carefully evaluated by means of 
two different approaches. On the one hand, the 
WD program provides the standard error associ- 
ated with each adjusted parameter. On the other 
hand, the error can also be estimated as the r.m.s. 
scatter of all the parameter sets corresponding to 
the iterations between the first solution and the 
last iteration (about 100). This is, in general, 
larger than the r.m.s. scatter of the solutions 
alone, i.e., those that fulfill the criterion previ- 
ously described. We conservatively adopted the 
uncertainty as twice whichever of the two error es- 
timates was the largest. The r.m.s. scatters of the 
residuals in the light curve fit are 0.011 (n — 264), 
0.006 (n = 303) and 0.005 mag (n = 303) for U, 
B, and V, respectively. Figure 1 shows the light 
curve fits to the observed U, B, and V differen- 
tial photometry and the corresponding residuals, 
where no systematic trends are observed. A 3D 
diagram of V380 Cyg (drawn to scale) is shown in 
Figure 3, where the tidal deformation of the larger 
primary component near periastron becomes evi- 
dent. As shown, the components are nearly spher- 
ical in shape at apastron. 

The results of the light curve fitting were com- 
bined with the velocity semi-amplitudes from 
the spectroscopic analysis (Paper I) and the 
UV/optical spectrophotometry study to obtain 
the absolute dimensions and radiative parameters 
listed in Table 4. The most recently published 
complete light curve study dates back to Hill & 



Batten (1984). These authors used the light 
program to analyze four light curves from various 
authors and with different qualities. In general, 
they determined a fractional radius for the pri- 
mary component that is significantly larger than 
our value, and for the secondary component this 
situation is reversed but the difference is much 
smaller. The reason for such a discrepancy prob- 
ably originates in the model used in their light 
curve solution, since the tidal distortion of the 
primary component is strong and also varies due 
to the eccentric orbit. In this sense, WD relics 
on a much more physically rigorous model than 
LIGHT and thus, our parameters ought to be pre- 
ferred. Moreover, the fractional radius of the pri- 
mary component is very well constrained from the 
shape of the out-of-eclipse curve near periastron 
(0?15), a region which is not well covered in the 
light curves analyzed by Hill & Batten (1984). 
Finally, these authors adopt an effective tempera- 
ture for the primary component (24 500 K), based 
on Johnson photometry which is clearly incon- 
sistent with our detailed analysis. However, their 
temperature ratio is in good agreement with our 
result. 

An inspection of the masses and radii obtained 
reveals a binary system composed of two mas- 
sive stars with rather different physical proper- 
ties: the more massive primary has evolved be- 
yond the main sequence whereas the lower-mass 
secondary star is still near to the zero-age main 
sequence (ZAMS). This situation makes V380 Cyg 
an extremely important astrophysical tool since it 
allows the stellar internal structure to be probed 
during a fast-evolving stage. Our determinations 
of masses and radii have accuracies of about 4 
percent and 2 percent, respectively. This is not 
as good Andersen (1991) suggests is required for 
a critical evaluation of stellar evolutionary mod- 
els. However, because of the short evolutionary 
timescales of the primary component of V380 Cyg 
and the strong dependence on the convection pa- 
rameters, this level of accuracy does allow for im- 
portant tests of the models, as we will show in §6. 

5. Times of Minima and Apsidal Motion 
Study 

V380 Cyg is a challenging system from the 
viewpoint of times of minima determinations: 
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both eclipses have long durations (more than 24 
hours), they have nearly- flat bottoms, and their 
shape is perturbed by the wave-like variation seen 
in the light curves. Therefore, the calculation of 
the exact time of central-eclipse is not straightfor- 
ward. Because complete coverage of one eclipse 
event is not possible during a single night, we de- 
cided to combine those observations correspond- 
ing to a similar epoch into a light curve and derive 
the times of the eclipses from comparison with the 
best available synthetic model. In our particular 
case, the most natural way of performing such op- 
eration is by considering the APT-PhocnixlO and 
the FCC 0.8-m APT data separately. The peri- 
ods span about 1 and 2.5 years and have around 
100 and 200 measurements per filter, respectively. 
An additional epoch was retrieved from obser- 
vations carried out during 1977/78 by Dorren 
& Guinan (unpublished) at Pahlavi Observatory 
(Shiraz, Iran). The exact time of minimum deter- 
mination was then done by fitting synthetic light 
curves computed with all the parameters listed 
in Table 3 to the actual observations at three 
different epochs. The resulting eclipse times are 
presented in Table 5. In addition, Table 5 lists 
the eclipse timings from Kron (1935), Semeniuk 
(1968), and Battistini et al. (1974). The mean 
HJD of each data set was adopted as the reference 
epoch. The secondary eclipse published by Seme- 
niuk (1968) was given a much lower weight in the 
calculations because of the poor coverage of the 
eclipses and thus the high uncertainty in the time 
determination. 

The eclipse timings listed in Table 5 allow 
the determination of linear ephemerides indepen- 
dently for primary and secondary eclipses with the 
following results (numbers between parentheses in- 
dicate errors in the last digits and E is the number 
of cycles): 

T Mi „i = HJD2441256. 544(6) + 12.425719(14) E 
T M inii = HJD2441261.625(5) + 12.425501(7) E 

The difference in the apparent periods is of course 
a clear indication of the presence of a significant 
apsidal motion. The complete determination of 
the apsidal motion rate of V380 Cyg was accom- 
plished in two steps. First, we computed the cor- 
responding time differences TMinii — ^Mini for each 
epoch and a linear least-squares fit to the data 



clearly showed the presence of a temporal decre- 
ment of 1.74 x 10~ 5 . The observations of Seme- 
niuk (1968) were not used because of the previ- 
ously mentioned lower accuracy. The result of the 
linear fit is shown in Figure 4, where the value of 
Semeniuk's measurement is plotted with an open 
circle. The measured slope corresponds to an ap- 
sidal motion rate of 0.0084 degrees per cycle, or a 
period of U = 1460 years, when an orbital eccen- 
tricity of 0.234, an inclination of 82.4 degrees, and 
a reference longitude of periastron of uo Q = 128 
degrees (representative of the average epoch) are 
adopted. In computing the apsidal motion rate, 
terms no higher than e 2 were used in the relevant 
equations. The slope determination appears to de- 
pend critically on Kron's (1935) early data point. 
However, when this measure is excluded, the re- 
sulting slope is still very similar (1.68 x 10~ 5 ) and 
within the expected uncertainties. 

As a second step, because of the relatively 
large orbital eccentricity, the method described by 
Gimenez & Garcia-Pelayo (1983), with equations 
revised by Gimenez & Bastero (1995), was used 
for a more accurate calculation of the apsidal mo- 
tion rate. When adopting the orbital inclination 
(82.4±0.2 deg) and the eccentricity (0.234±0.008) 
as derived from the light curve analysis (see Ta- 
ble 3), we obtained an apsidal motion rate of 
0.0083 ± 0.0006 degrees per cycle, equivalent to 
cj bs = 24.0 ± 1.8 degrees per 100 years or an apsi- 
dal motion period of U = 1490 ± 120 years. These 
values are in excellent agreement with previous 
determinations and the position of the periastron 
predicted with these elements for the epoch of the 
radial velocity curve of Paper I, HJD 2449500, is 
uo = 132. 4 ±1.5 degrees also in excellent agreement 
with that analysis (w = 133.2 ± 3.2). 

6. Evolutionary Model Study 

Stellar evolutionary theory provides a descrip- 
tion of the interior structure and the observable 
properties of a star of given initial mass and chem- 
ical composition as a function of age. Signifi- 
cant progress has been made in recent years in 
understanding the physical processes that govern 
the structure and evolution of stars, including the 
computation of accurate nuclear reaction rates, 
neutrino emission rates, and atmospheric opaci- 
ties. 
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Modeling of the convective energy transport 
mechanism is probably one of the major deficien- 
cies of the current stellar structure and evolution 
models (see Gimenez, Guinan, & Montesinos 1999 
and references therein). Because of the extreme 
mathematical complexity of an accurate treat- 
ment, convection is still described through a very 
simplified and phenomcnological theory: the mix- 
ing length theory (MLT). When comparing the 
predictions of the models for stars with convective 
nuclei, an extra amount of convection (overshoot- 
ing) is needed to reproduce the observed proper- 
ties. The overshooting parameter (a ov ), i.e., the 
size in pressure-scale height of the over-extension 
of the convective nucleus, cannot be theoretically 
computed and thus needs to be determined from 
observations. The significance of overshooting is 
probably one of the largest uncertainties in MLT, 
as well as the possible dependence of this param- 
eter on mass, evolution or chemical composition. 
Some attempts (Canuto & Mazzitelli 1991, 1992; 
Canuto 1999) have been directed towards the de- 
velopment of an essentially parameter-free convec- 
tion model, but the MLT is still being widely used 
with fairly good results in most applications. 

The eclipsing binaries are powerful tools for 
testing stellar structure and evolution models, 
since the fundamental properties of the compo- 
nents (masses, radii, luminosities, etc.) can be 
accurately determined from observations. Several 
studies testing the significance of the convective 
overshooting in the stellar nucleus using eclips- 
ing binary data (Andersen, Nordstrom, & Clausen 
1990; Pols et al. 1997; Ribas et al. 1999) have 
been published. Although these studies clearly 
show an improvement when overshooting is in- 
cluded, more quantitative conclusions are difficult 
to access. This is primarily because the effects 
of overshooting are only observationally relevant 
for evolved stages, especially beyond the TAMS. 
However, no evolved early-type detached eclipsing 
binaries with accurate determination of absolute 
dimensions were available to date because of ob- 
servational selection effects. V380 Cyg is in this 
sense an important and unique system. 

Evolutionary tracks (specifically generated for 
V380 Cyg) for three different overshooting param- 
eters (a v = 0.2, 0.4, and 0.6) were kindly pro- 
vided by A. Claret, for the comparison of the ob- 
served stellar properties with the theoretical pre- 



dictions. The input physics for these models is 
identical (except for the enhanced mixing) to that 
described in Claret (1995), i.e., OPAL opacities, 
mass loss beyond 10 Mq, mixing length parame- 
ter a p — 1.52, modern equations of state, exten- 
sive nuclear network, etc. The initial masses of 
the evolutionary tracks were selected so that they 
would agree with the observed masses of the com- 
ponents at their current evolution stage. A metal 
abundance of Z = 0.012 was adopted as discussed 
in §3. The helium abundance (Y) that yields the 
best match to the observed properties of the es- 
sentially unevolved secondary component (not af- 
fected by the value of the overshooting parame- 
ter, see below) was found to be Y = 0.26 ± 0.02. 
This value is in good agreement with that expected 
from the metal abundance of V380 Cyg and an av- 
erage helium-to-metal enrichment law (see Ribas 
et al. 2000a). A logT e ff — log g plot showing the 
system components and the tracks computed for 
the adopted chemical composition is presented in 
Figure 5. 

The effect of adopting different values of the 
overshooting parameter in the early post-ZAMS 
evolutionary stage of the secondary component is 
very small. Therefore, its observed properties can 
be well reproduced by models with wide ranges 
of convective overshooting. The evolved primary 
component, however, with logg = 3.15, presents 
a quite different situation. This star provides an 
excellent test of convective overshooting. Indeed, 
the evolutionary track with a ov — 0.2 predicts an 
effective temperature for the surface gravity of the 
primary component that is ~ 3000 K (~ 7a) lower 
than the observed value. Even with a ov = 0.4, the 
star is predicted to have an effective temperature 
2000 K lower than observed. A value of a ov = 0.6 
is required to achieve consistency between model 
predictions and observations, as shown in Figure 5. 
In this case, the primary component is predicted 
to be close to the phase of hydrogen shell ignition. 

Although not explicitly shown in Figure 5 for 
the sake of clarity, evolutionary tracks for the ob- 
served mass plus and minus 1 a were also available. 
The M + la model for the primary component 
yields effective temperatures at a given log g only 
about 0.009 dex (w 500 K) larger than the model 
for the observed mass. Therefore, even when con- 
sidering the error bars, the location of the primary 
component is not compatible with evolutionary 
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tracks including mild overshooting. A range of 
overshooting parameters within a ov « 0.6 ± 0.1 is 
possible when all the observational errors (mass, 
log g, chemical composition, temperature ratio, 
etc) are taken into account in the analysis. 

Even though the effective temperature is a well- 
determined quantity in our analysis, it might still 
be affected by systematics in the model atmo- 
sphere fluxes or in the IUE observations. How- 
ever, our result favoring a large amount of convec- 
tive core overshooting is immune to any of these 
effects, because the temperature ratio of the com- 
ponents is accurately determined from the light 
curve analysis. The good agreement obtained for 
the secondary component excludes the possibility 
of resorting to a different value of Y or an er- 
ror in the zero-point of the temperature scale for 
explaining the temperature difference. The con- 
clusion of the analysis is based on the tightly- 
constrained relative location of the components 
in the \ogg — logT ff diagram (i.e., effective tem- 
perature ratio) rather than on the absolute value 
adopted for the effective temperatures. 

It is important to note that in a detached sys- 
tem such as V380 Cyg the ages predicted by the 
models for the two components should agree. This 
restriction is not very critical for V380 Cyg, how- 
ever, because the secondary component is close to 
the ZAMS and, since the observational errors in 
both the mass and the radius is of the order of 
2-3 percent, a wide age range is compatible with 
the stellar parameters. This is shown in Table 6, 
where the model predictions (and the errors) are 
given for the three overshooting values considered. 
The age, effective temperature, and surface grav- 
ity for the primary component are presented in the 
table. Also included is the age of the secondary 
component. These parameters were interpolated 
in the computed tracks at the observed values of 
log g for both components. The evolutionary ages 
for the primary and secondary components are 
identical within the errors and no significant age- 
dependent variations are apparent for the various 
values of the overshooting parameter. Therefore, 
the coeval criterion for the stars does not favor any 
particular overshooting parameter but it confirms 
the good agreement between observed and model 
predicted physical properties. The age of the sys- 
tem was found to be r = 25.5 ± 1.5 Myr when the 
best-fitting a ov « 0.6 was adopted. 



A detailed analysis was also carried out with 
evolutionary tracks of solar metal content (Z = 
0.02, Y = 0.28). From the locations of the stars 
in the log g — logT e ff we found that a ov w 0.6 
still clearly yields the best fit, although the com- 
ponents are not compatible with coeval evolution 
(a small age differential of ~ 3 Myr between the 
components becomes evident). 

It has been suggested that the effects of ro- 
tation on the structure and evolution of massive 
stars could mimic those of overshooting. We used 
evolutionary tracks for V380 Cyg provided by A. 
Claret (the input physics and approximations are 
described in detail in Claret 1999) to explore this 
possibility. Models that include rotation do not 
show significant differences in the logT e ff — \ogg 
diagram when compared to the non-rotating mod- 
els, especially in the evolved stages (a similar con- 
clusion was reached by Deupree 1998). Thus, it 
appears that this effect cannot mitigate the rel- 
atively high value of the overshooting parameter 
found here. 

In practice, the inclusion of rotation in the mod- 
els only has significant effect when the rotational 
velocity of the star is fast and approaches the 
critical value (break up). However, V380 Cyg 
is a slow rotator (P « 12 days). The initially 
faster-rotating eclipsing binary components are 
braked by tidal effects until they become synchro- 
nized with the orbital period. Because the lat- 
ter is usually quite long for detached massive or 
intermediate-mass systems, the components spin 
more slowly than single stars of the same mass. 
Therefore, evolutionary models that do not in- 
clude rotation can more confidently be used in 
eclipsing binaries with massive components than 
with more rapidly rotating single stars (in clus- 
ters, for example). 

The determination of the apsidal motion rate 
of V380 Cyg permits an independent test of stel- 
lar model predictions. The apsidal motion rate of 
an eccentric binary system can be computed from 
theory as the sum of the contribution of general 
relativity (GR) and classical (CL) effects (see e.g., 
Guinan & Maloney 1985). In a system such as 
V380 Cyg, with a distorted component, the CL 
contribution to the apsidal motion is expected to 
dominate the GR term. The classical apsidal mo- 
tion rate is a function of the orbital (P, e) and 
physical (Mp ; s, rp t s, w ro tp,s) properties of the sys- 



9 



tern, all determined from the analysis of the light 
and radial velocity curves. But in addition, the 
classical apsidal motion depends on the internal 
mass distribution of the stars. A star's internal 
mass distribution is parameterized in polytropic 
models by the internal structure constant k 2 - In 
this theory, k 2 is proportional to the ratio of the 
mean and the central densities of the star. The 
values of k 2 can be computed from stellar evolu- 
tion codes (see Claret & Gimenez 1992). 

An observed apsidal motion rate of w(obs) = 
24.0 ± 1.8 °/100 yr was obtained in §5. For com- 
parison with the model predictions, k 2 can be com- 
puted from w(obs) (after correcting for the gen- 
eral relativistic effect, uqr — 2 °/100 yr) by using 
the expressions in Claret & Gimenez (1993) and 
by assuming pseudo-synchronization. Addition- 
ally, several studies (e.g., Stothers 1974; Claret 
& Gimenez 1993) indicate that a small correction 
may be needed to include the effects of stellar ro- 
tation in the calculation of k 2 . 

Taking these effects into consideration, we com- 
puted a systemic mean observational value of 
logfc 2 (obs) = -2.89 ± 0.04. The contribution 
of the secondary component to the mean inter- 
nal structure constant is negligible (less than 1 
percent) because of its small fractional radius. 
Thus, we assume log/c 2 p(obs) = —2.89 ± 0.05. 
The specially constructed evolutionary models for 
V380 Cyg, kindly provided by A. Claret, yield the 
values listed in Table 6. As can be seen, a value 
of a ov as large as 0.6 ± 0.1 is clearly needed to fit 
the observed apsidal motion rate of the system. 

It is remarkable that the overshooting param- 
eter of a ov « 0.6 ±0.1 is supported by two inde- 
pendent approaches, namely the investigation of 
the apsidal motion and the location of the compo- 
nents in the theoretical H-R diagram. The value 
indicated is significantly larger than that currently 
adopted by most theoretical models. Some brief 
comments on the physical implications of this re- 
sult are left for §8. In addition, a more exten- 
sive study on stellar core convection using all the 
available data from eclipsing binaries is the subject 
of a forthcoming paper (Ribas, Jordi, & Gimenez 
2000b). 



7. Tidal Evolution 

V380 Cyg has an eccentric orbit and the com- 
ponent stars appear to have rotational velocities 
synchronized with the orbital velocity at perias- 
tron, the so called pseudo-synchronization (see Ta- 
ble 4). Pseudo-synchronization is commonly found 
in eccentric binary systems because in most cases 
the circularization timescale is much larger than 
the synchronization timescale (Claret, Gimenez, & 
Cunha 1995). To compare the observations with 
the theoretical predictions we employed the tidal 
evolution formalism of Tassoul (1987, 1988). The 
expressions used, which were taken from Claret et 
al. (1995), relate the time variation of the eccen- 
tricity (e) and the angular rotation of the com- 
ponents (f2) with the orbital and stellar physical 
properties (mass ratio, orbital period, mass, lu- 
minosity, radius and gyration radius). It is im- 
portant to integrate the differential equation that 
governs the variation of the eccentricity and the 
angular velocity with time rather than use a sim- 
ple timescale calculation. This is especially advis- 
able for evolved systems like V380 Cyg because 
the components may have undergone large radius 
and luminosity changes that have an important 
effect over the circularization and synchronization 
timescales. 

Thus, we integrated the differential equations 
along evolutionary tracks adopting the observed 
masses of the V380 Cyg components. In these cal- 
culations we assumed a constant orbital period. 
As boundary conditions, we assumed the eccen- 
tricity to be the observed value (e = 0.234) at the 
present age of the system (r = 25.5 Myr), and 
the angular rotation of the stars to be normal- 
ized to unity at the ZAMS. The calculations were 
made with the evolutionary models at Z — 0.012, 
Y = 0.26, and a ov = 0.6. The predictions of the- 
ory agree well with a system that has reached syn- 
chronism for both components but that has not yet 
circularized its orbit. These calculations indicate 
that the primary component reduced the differ- 
ence between its angular velocity and the pseudo- 
synchronization angular velocity to ~ 0.1 percent 
of the ZAMS value at the early age of 10 Myr, 
and the same occurred to the secondary compo- 
nent at an age of ~ 20 Myr. 

From the boundary condition it is inferred that 
the initial (ZAMS) orbital eccentricity was e « 
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0.3. A comparison of this initial value with the 
current observed one indicates a net decrease of 
Ae w 0.07 in about 25 Myr. The theory also pre- 
dicts that circularization should occur on a rela- 
tively short timescale of w 1.5 x 10 4 yr when the 
radius of the primary star rapidly increases shortly 
before the initiation of core helium burning. The 
instantaneous (i.e., with the current parameters) 
circularization timescale is « 10 8 yr, much larger 
than the value computed from the integration of 
the differential equation because the primary com- 
ponent is located in a very fast evolving stage. 
Although tidal evolution theory predicts orbital 
circularization for V380 Cyg in 1.5 x 10 4 yr, the 
Roche lobe constraints are not included in the 
modeling. If these are included, the primary star 
should fill its Roche lobe in about 10 4 yr. This re- 
sults in Roche lobe overflow and the initiation of 
mass transfer. Once V380 Cyg reaches the semi- 
detached stage, the binary should circularize very 
rapidly due to strong tidal interactions and fric- 
tion. 

8. Conclusions 

This study complements the spectroscopic in- 
vestigation of V380 Cyg published in Paper I. 
New differential photometry has been collected in 
two different epochs, resulting in about 300 ob- 
servations in each of the UBV bandpasses. The 
photometry has been modeled using an improved 
version of the Wilson-Devinney program that in- 
cludes recent model atmospheres. The best fit to 
the light curves was achieved for the stellar and 
orbital properties listed in Table 3. The effective 
temperatures of the stars were determined using 
UV spectrophotometry and optical photometry. 
Also, four new eclipse timings were added to the 
published ones and a complete re-evaluation of the 
apsidal motion rate of V380 Cyg was carried out. 

The determined physical properties (Table 4) 
of the V380 Cyg components were compared with 
evolutionary models (with the same input physics 
as in Claret 1995) computed for the observed 
masses, and with different values of convective 
overshoot (a ov = 0.2, 0.4, and 0.6). The fractional 
abundances of metals and helium were estimated 
to be Z = 0.012 ± 0.003 (from UV/optical spec- 
trophotometry) and Y = 0.26 ± 0.02 (from evo- 
lutionary model fit to the secondary component) 



respectively. The best agreement was found for 
an overshooting parameter of a ov « 0.6, for which 
the primary component is predicted to be located 
near the blue point of the MS hook. The excel- 
lent agreement obtained for the secondary com- 
ponent excludes the possibility of resorting to a 
different Y value or an error in the temperature 
scale zero point for explaining the location of the 
primary component, because the temperature ra- 
tio of the components is well-constrained from the 
light curve analysis. The apsidal motion study of 
V380 Cyg also indicates that a value of a ov w 0.6 
is necessary to fit the observed apsidal motion rate 
of the system. The study of the locations of the 
components in the \ogg — logT e g diagram and 
the internal structure (from apsidal motion) lead 
to mutually consistent results. However, we cau- 
tion that the relatively large value of overshoot- 
ing should be restricted to the mass, evolutionary 
stage, and chemical composition of the primary 
component of V380 Cyg. 

Our convective overshooting result has a more 
profound implication than merely providing a bet- 
ter fit to the V380 Cyg data. We have shown that 
the observed physical properties of V380 Cyg can- 
not be described by evolutionary models with the 
most commonly used physical ingredients (includ- 
ing mild overshooting). The physical difference 
between these standard models and our better- 
fitting models (with enhanced convective over- 
shooting) is that the latter have larger convective 
cores. Thus, our more general conclusion is that 
massive stars have larger convective cores and, 
therefore, are more centrally condensed than pre- 
dicted by standard theory. Note that convective 
overshooting is not the only physical process that 
might produce such a result. Further increases in 
the opacity (as it was demonstrated when substi- 
tuting the LAOL by the newer OPAL opacity ta- 
bles), diffusive mixing, turbulence or other mech- 
anisms may certainly lead to stellar models with 
larger and denser convective cores, relaxing the 
need for our large overshooting parameter. In the 
meantime, with the current physical ingredients, 
our results indicate that models should consider 
a moderate amount of convective overshooting for 
matching the observed properties of real massive 
stars. 

An investigation of the tidal evolution of 
V380 Cyg was also carried out. The numeri- 



11 



cal integrations of the relevant differential equa- 
tions indicate a binary system with components 
in an eccentric orbit that have reached pseudo- 
synchronism. From the theory of Tassoul (1987, 
1988), it appears that the orbit circularization 
should occur on a timescale of 10 4 years. This 
is on the same timescale as it takes the primary 
component to fill its Roche lobe and start mass 
transfer. 

This paper clearly demonstrates the impor- 
tance of specific eclipsing binaries as "astrophysi- 
cal laboratories" . In the present case of V380 Cyg, 
fundamental problems of stellar structure and evo- 
lution are addressed. Early-type massive systems 
with at least one evolved component constitute a 
very important source of observational data for 
testing stellar structure and core convection. Such 
systems are scarce (due to strong observational se- 
lection effects) and gathering the required photo- 
metric and spectroscopic data is not easy (because 
of the typically long orbital periods - 10 days or 
more). The relative faintness of the secondary 
component also presents difficulties in the spec- 
troscopic aspects of the problem. In spite of the 
observational challenges, these systems (although 
rare) have proved to be crucial tools and their 
study should be a priority in stellar astrophysics. 

Some examples of early-type eclipsing bina- 
ries suitable for testing core convection are: 
V1756 Cyg, V453 Cyg and V346 Cen (see Ribas 
et al. 1999). V1765 Cyg has a very massive 
and evolved component, already in the supergiant 
stage. V453 Cyg and V346 Cen have primary com- 
ponents of similar mass to that of V380 Cyg but 
appear to be at a less evolved stage (log<? « 3.8). 
Other suitable binaries for testing models are 
those belonging to the £-Aur class of eclipsing 
binaries (Schroder, Pols, & Eggleton 1997). These 
are systems with at least one evolved member in 
the core helium burning or in the supergiant stage. 
The long periods (hundreds or even thousands of 
days) make the observations challenging. Hope- 
fully, there will be many more potentially interest- 
ing systems to study. The long term monitoring 
programs associated with (future) space missions 
and the ground-based systematic surveys (e.g., the 
Sloan Digital Survey and the OGLE experiment 
towards the Galactic Bulge) are expected to yield 
valuable results, both from the detection point of 
view and from the densely-covered light curves 



that they will provide. Finally, extensive photom- 
etry of extragalactic eclipsing binaries (in LMC, 
SMC, and M31) are becoming available as a result 
of several microlensing (Grison et al. 1995; Alcock 
et al. 1997; Udalski et al. 1998) and variable star 
(Kaluzny et al. 1998, 1999; Stanek et al. 1998, 
1999) surveys. Such a wealth of data guarantees 
a nearly inexhaustible source of eclipsing systems 
meeting the demanding requirements for critical 
model analysis, not only for solar-type abundances 
but also for other chemical compositions. 
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Daniel M. Popper who died last September before 
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Fig. 1. — Light curve fit to the observed U, B and V Comparison— Variable (Cp— V) differential photometry 
of V380 Cyg. Also shown are the Observed— Computed (O— C) residuals. The orbital phase was computed 
according to the following ephemeris: T M ini = HJD2441256.544 + 12.425719 E. 



15 



-9.6 




-1 1.0 t i i i i □ 

1100 2000 3000 4000 5000 6000 

Wavelength (A) 



Fig. 2.— Kurucz atmosphere model fit (solid line) to IUE SWP+LWP spectrum (dots) of V380 Cyg. The 
cross symbols indicate wavelength regions that were not used in the fits chiefly because of the presence 
of contaminating interstellar medium features. The filled circles and triangles are Johnson and Stromgren 
fluxes respectively. The residuals of the fit are shown in the lower part of the figure. 
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Fig. 3. — 3-D picture of V380 Cyg at periastron and apastron. Notice the slight deformation of the primary 
component during periastron. 
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Fig. 4. — Apsidal motion calculation through a linear fit to observed time difference between primary and 
secondary minima. The slope of the straight line is 1.74 x 10~ 5 . The individual points have been computed 
from the times of minima listed in Table 5. Semeniuk's (1968) value is plotted as an open circle and was not 
used in the fit (see text). 
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Fig. 5. — log g — logTcff plot of V380 Cyg. Evolutionary tracks for the primary (P) and secondary (S) 
components computed with a ov = 0.2, 0.4, and 0.6 are shown. 
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Table 1 

Differential photometry for V380 Cyg in the UBV bandpasses. 
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8901 





.263 


49891 


8915 





.273 


49891 


.8925 





269 


49891 


8939 


0.277 


49891 


.8954 





.272 


49892 


8896 





286 


49892 


8910 





295 


49892 


.8920 





.290 


49892 


.8934 





294 


49892 


.8949 


0.291 


49893 


.8815 





.231 


49893 


8829 





242 


49893 


8838 





238 


49893 


.8853 





.237 


49893 


.8868 





236 


49895 


8369 


0.250 


49895 


.8379 





.243 


49895 


8393 





249 


49895 


8408 





212 


49896 


.8945 





.298 


49896 


.8960 





.300 


49896 


8970 


0.303 


49896 


.8984 





.302 


49896 


8999 





299 


49898 


8476 





248 


49898 


.8490 





.255 


49898 


.8499 





255 


49898 


8514 


0.254 


49898 


.8529 





.252 


49899 


8734 





246 


49899 


8748 





263 


49899 


.8758 





.255 


49899 


.8772 





260 


49899 


8787 


0.255 


49900 


.8428 





.243 


49900 


8443 





255 


49900 


8467 





260 


49901 


.8803 





.239 


49901 


.8817 





.247 


49901 


8828 


0.247 
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Table 1 — Continued 





HJD- 


V-C 


HJD- 


V 


-c 


HJD- 


V 


-c 


HJD- 


V-C 


HJD- 


V-C 


HJD- 


V 


-C 


2400000 






2400000 






2400000 






2400000 






2400000 






2400000 






49901 


8842 





249 


49901 


8857 





240 


50237.9721 





242 


50237.9735 


0.250 


50237 


9750 





247 


50251 


9588 





258 


50252 


9619 





351 


50252 


9649 





348 


50253.9529 





239 


50253.9543 


0.246 


50253 


9558 





257 


50254 


9481 





228 


50254 


9496 





235 


50254 


9511 





223 


50259.9446 





248 


50259.9460 


0.258 


50259 


9475 





250 


50396 


5739 





247 


50396 


5753 





253 


50396 


5768 





252 


50713.6780 





252 


50713 


6794 





257 


50714 


6783 





228 


50714 


6798 





236 


50718 


7015 





249 


50718 


7029 





259 


50718 


7044 





248 


50721 


6658 





239 


50721 


6672 





254 


50721 


6688 





243 


50725 


6386 





266 


50725 


6401 





275 


50725.6416 





268 


50726.6661 





223 


50726 


6676 





232 


50726 


6691 





229 


50727.6719 





238 


50727 


6734 





246 


50727.6749 





240 


50728.6342 





256 


50728 


6356 





263 


50728 


6372 





259 


50730.6658 





256 


50752 


5807 





237 


50752 


5821 





211 


50752 


5836 





235 


50753 


5973 





257 


50753 


5987 





259 


50753 


6002 





250 


50754 


5966 





334 


50754 


5997 


0.336 


50755 


5958 





258 


50755 


5973 





264 


50755 


5988 





258 


50756 


5952 





252 


50756 


5967 





266 


50756 


5982 





257 


























V 


47415 


8013 





217 


47417 


7500 





204 


47418 


7531 





227 


47420 


7519 





221 


47421 


7429 





196 


47422 


7256 





215 


47423 


7161 





229 


47424 


7076 





297 


47427 


7040 





207 


47428 


7033 





208 


47429 


7106 





208 


47431 


7010 





256 


47432 


7083 





239 


47433 


7389 





186 


47436 


7233 





268 


47441 


7282 





209 


47443 


7104 





219 


47459 


6380 





204 


47461 


6359 





287 


47466 


6301 





207 


47469 


6220 





311 


47470 


6211 





200 


47472 


6245 





217 


47475 


6085 





218 


47476 


6031 





221 


47479 


5852 





205 


47482 


5705 





228 


47483 


5706 





194 


47485 


5705 





221 


47486 


5703 





288 


47493 


5622 





234 


47494 


5554 





288 


47495 


5612 





201 


47497 


5610 





222 


47498 


5546 





227 


47500 


5611 





216 


47508 


5560 





193 


47618 


0173 





247 


47619 


0159 





246 


47620 


0132 





189 


47622 


0135 





225 


47623 


0106 





256 


47624 


0024 





231 


47626 


9985 





202 


47629 


9922 





218 


47631 


9828 





215 


47632 


9770 





187 


47635 


9776 





288 


47636 


9714 





216 


47641 


9575 





218 


47646 


9470 





226 


47647 


9447 





264 


47648 


9401 





219 


47649 


9361 





223 


47651 


9321 





213 


47664 


8964 





212 


47665 


8909 





204 


47668 


8804 





234 


47669 


9391 





194 


47670 


8765 





197 


47671 


9653 





225 


47672 


8840 





280 


47673 


8940 





225 


47677 


8624 





203 


47680 


8474 





301 


47683 


8784 





206 


47686 


8887 





216 


47687 


8811 





208 


47688 


8657 





199 


47689 


8647 





215 


47690 


8620 





199 


47691 


8520 





216 


47693 


8433 





222 


47694 


8597 





191 


47695 


8595 





208 


47699 


8585 





214 


47701 


8618 





207 


47779 


7548 





293 


47780 


7559 





232 


47788 


7323 





199 


47789 


7426 





211 


47793 


7278 





208 


47794 


7533 





200 


47797 


7467 





291 


47799 


7394 





210 


47800 


7425 





209 


47812 


6939 





206 


47813 


6851 





201 


47814 


6905 





211 


47816 


6890 





216 


47817 


6812 





253 


47818 


6881 





198 


47826 


6739 





207 


47827 


6683 





203 


47833 


6399 





221 


47834 


6378 





293 


47835 


6322 





214 


47836 


6230 





219 


47838 


6264 





213 


47840 


6143 





209 


49864 


9461 





200 


49864 


9475 





205 


49864 


9485 





202 


49864 


9499 





208 


49864 


9514 





203 


49866 


9332 





226 


49866 


9346 





219 


49866 


9356 





224 


49866 


9370 





228 


49866 


9385 





221 


49867 


9054 





282 


49867 


9069 





280 


49867 


9079 





279 


49867 


9093 





280 


49867 


9108 





279 


49868 


9410 





201 


49868 


9425 





207 


49868 


9434 





207 


49868 


9448 





210 


49868 


9463 





199 


49869 


9476 





194 


49869 


9490 





201 


49869 


9500 





196 


49869 


9514 





200 


49869 


9529 





197 


49870 


9379 





213 


49870 


9393 





213 


49870 


9403 





214 


49870 


9417 





215 


49870 


9432 





213 


49871 


9392 





224 


49871 


9406 





229 


49871 


9416 





222 


49871 


9430 





226 


49871 


9445 





226 


49872 


9385 





246 


49872 


9399 





256 


49872 


9409 





254 


49872 


9423 





253 


49872 


9438 





253 


49873 


9201 





217 


49873 


9216 





226 


49873 


9226 





211 


49873 


9240 





221 


49873 


9255 





217 


49876 


8986 





201 


49876 


9000 





201 


49876 


9010 





201 


49876 


9024 





205 


49876 


9039 





195 


49878 


9427 





216 


49878 


9441 





221 


49878 


9451 





212 


49878 


9465 





212 


49878 


9480 





216 


49879 


9266 





299 


49879 


9280 





303 


49879 


9290 





299 


49879 


9304 





307 


49879 


9319 





302 


49880 


9016 





214 


49880 


9031 





222 


49880 


9040 





219 


49880 


9055 





221 


49880 


9069 





228 


49883 


9258 





209 


49883 


9272 





219 


49883 


9282 





218 


49883 


9296 





224 


49883 


9311 





212 


49884 


8756 





283 


49884 


8770 





290 


49884 


8780 





286 


49884 


8794 





295 


49884 


8809 





285 


49885 


9684 





208 


49885 


9699 





217 


49885 


9708 





215 


49885 


9722 





222 


49885 


9737 





217 


49886 


9113 





210 


49886 


9127 





212 


49886 


9137 





209 


49886 


9151 





215 


49886 


9166 





208 


49887 


8954 





206 


49887 


8968 





210 


49887 


8978 





204 


49887 


8993 





209 


49887 


9007 





203 


49888 


9282 





199 


49888 


9296 





203 


49888 


9306 





205 


49888 


9320 





201 


49888 


9335 





199 


49889 


9535 





198 


49889 


9549 





206 


49889 


9559 





202 


49889 


9573 





201 


49889 


9588 





200 


49890 


8934 





209 


49890 


8948 





216 


49890 


8958 





207 


49890 


8972 





203 


49890 


8987 





201 


49891 


8902 





222 


49891 


8917 





228 


49891 


8927 





221 


49891 


8941 





228 


49891 


8956 





224 


49892 


8897 





241 


49892 


8912 





245 


49892 


8921 





246 


49892 


8935 





251 


49892 


8950 





247 


49893 


8816 





193 


49893 


8830 





195 


49893 


8840 





193 


49893 


8854 





200 


49893 


8869 





195 


49895 


8356 





206 


49895 


8370 





207 


49895 


8380 





203 


49895 


8394 





204 


49895 


8409 





204 


49896 


8946 





253 


49896 


8961 





260 


49896 


8971 





254 


49896 


8986 





260 


49896 


9001 





252 


49898 


8477 





204 


49898 


8491 





209 


49898 


8501 





207 


49898 


8515 





208 


49898 


8530 





207 


49899 


8736 





202 


49899 


8750 





213 


49899 


8760 





209 


49899 


8774 





211 


49899 


8789 





211 


49900 


8430 





208 


49900 


8444 





210 


49900 


8468 





211 


49901 


8804 





196 


49901 


8819 





205 


49901 


8829 





201 


49901 


8843 





206 


49901 


8858 





201 


50237 


9722 





205 


50237 


9737 





207 


50237 


9752 





208 


50251 


9590 





213 


50252 


9621 





306 


50252 


9635 





315 


50252 


9650 





307 


50253 


9530 





200 


50253 


9544 





213 


50253 


9559 





205 


50254 


9483 





184 


50254 


9497 





192 


50259 


9447 





204 


50259 


9462 





215 


50259 


9476 





209 


50396 


5740 





209 


50396 


5754 





209 


50396 


5769 





209 


50713 


6796 





212 


50714 


6785 





185 


50714 


6799 





193 


50718 


7016 





204 


50718 


7031 





214 


50718 


7046 





206 


50721 


6659 





199 


50721 


6674 





205 


50721 


6689 





205 


50725 


6388 





219 


50725 


6402 





232 


50725 


6417 





224 


50726 


6663 





181 


50726 


6677 





188 


50726 


6693 





185 


50727.6721 





202 


50727 


6735 





209 


50727 


6751 





196 


50728 


6343 





215 


50728 


6358 





215 


50728 


6373 





213 


50730 


6644 





210 


50730 


6659 





212 


50730 


6674 





209 


50752 


5823 





200 


50752 


5838 





194 


50753 


5974 





207 


50753 


5989 





220 


50753 


6004 





214 


50754 


5968 





294 


50754 


5998 





293 


50755 


5960 





213 


50755 


5974 





219 


50755.5990 


0.216 


50756 


5954 





211 


50756 


5968 





221 


50756 


5984 





217 
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Table 2 

Photometric properties and effective temperature determination of V380 Cyg 



Johnson a 


Str6mgren b 




Vilnius 


V = 5.68 ±0.02 


b -y = 0.031 


U - 


- P = 0.20 


B — V = -0.06 ±0.01 


mi = 0.039 


P - 


-Y = 0.37 


U -B = -0.76 ±0.02 


ci = 0.111 


Y - 


-V = 0.22 




(3 = 2.587 






E(B -V) = 0.18 
Toff = 22 700 K (F96) 
T c ff = 24 600 K (P80) 
Toff = 23 100 K (B81) 


E(B -V) = 0.19 
Toff = 20 200 K 


E(B 

T e ff 


- V) = 0.23 
= 20 900 K 



a Photometry from this work. E{B — V) computed using the 
intrinsic (U — B) vs (B — V) relation from Schmidt-Kalcr 1982 
(luminosity class III) and the standard reddening slope of E(U — 
B)/E(B — V) = 0.72. Effective temperatures obtained from T e g 
vs. (B-V)o calibrations of Flower 1996 (F96), Popper 1980 (P80) 
and B6hm-Vitense 1981 (B81). 

b Photomctry from Hauck & Mcrmilliod 1998. E(B - V) deter- 
mined from Crawford 1978 — as explained in Figueras, Torra, & 
Jordi 1991 and Jordi et al. 1997 — and temperature computed 
by means of the photometric grids of R. Napiwotzki 1998 (private 
communication), based on Kurucz ATLAS9 atmosphere models. 
Coefficients for the color excesses taken from Straizys 1992 and 
Crawford & Mandwewala 1976. 

c Photometry from Sudzius et al. 1992. Calibrations in Straizys 
(1992) used for the calculation of the color excess and the effective 
temperature. 
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Table 3 

Results from the light curve, radial velocity curve and UV/optical spectrophotometry 

analyses. 



Lp 
Ls 
Lp 
is 
Lp 
Ls 



Parameter Value 

Light curve 

Period (days) 12.425719 ± 0.000014 

Eccentricity 0.234 ± 0.006 

Inclination (deg) 82.4 ± 0.02 

lo (1989.0) (deg) 132.7 ±0.3 

lu (deg/100 yr) 24.0 ±1.8 

^ 0.96 ±0.02 

19.14 ±0.11 

u 

17.48 ±0.08 

B 

17.31 ±0.07 

r P 0.2478 ±0.0014 

rs 0.0631 ± 0.0008 

Radial velocity curve a 
K P (km s- 1 ) 94.5 ±1.5 

K s (km s- 1 ) 151.1 ±3.0 

g=H 0.626 ±0.041 

7 (kmV 1 ) 2.4 ±1.0 

a (R ) 59.2 ±0.8 

UV/Optical Spectrophotometry 
T cffP (K) 21 350 ± 400 

«microp (km s _1 ) 12 ± 1 

^microS 

(km s^ 1 ) 

[m/H] (dex) -0.44 ±0.07 

- V) (mag) 0.17 ±0.02 

log(i?p/d 2 ) -18.961 ±0.031 

a From Paper I. 
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Table 4 

Properties of the V380 Cyg components. 





Primary 


Secondary 


Spectral Type a 


B1.5 II-III 


B2 V 


Mass (M j 


11.1 ±0.5 


6.95 ±0.25 


Radius (R©) 


14.7 ±0.2 


3.74 ±0.07 


log g (cgs) 


3.148 ±0.023 


4.133 ±0.023 


Toff (K) 


21 350 ±400 


20 500 ±500 


log(L/L ) b 


4.60 ±0.03 


3.35 ±0.04 


M bol c (mag) 


-6.75 ±0.07 


-3.62 ±0.10 


V d (mag) 


5.74 ±0.02 


8.84 ±0.02 


d°(pc) 


1000 ±40 


M v f (mag) 


-4.79 ±0.10 


-1.69 ±0.10 


w rot s (km s" 1 ) 


98 ±4 


32 ±6 


w rotp— sync (km S ) 


99 ±2 


25 ±2 


p (g cm" 3 ) 


0.0049 ± 0.0003 


0.187 ±0.013 


pj (g cm" 3 ) 


3.09 ±0.24 


19.6 ±1.6 



a Hill & Batten 1984. 
b From L = AitR 2 <jT^ s . 

c From \og(L/ Lq) and adopting Mboi© = 4.75. 

d Computed from the out-of-eclipse magnitude and 

(£p/£s)lv- 

e Using \og(Rp/d 2 ) from the UV/optical spectropho- 
tometry fit and the observed Rp. 

'From M v = V - 3.1E(B - V) + 5 - 51ogd. Note 
that our values are completely consistent with what is 
expected for the spectral classification of the components. 
In addition, M\, \ — M v for both components is in excellent 
agreement with the bolomctric correction calibration of 
Flower (1996). 

g Computed from Lyubimkov et al. 1996. 

h Pseudo-synchronization rotational velocity. 

'Following Kopal 1959 and using the internal concen- 
tration parameters k 2 . 
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Table 5 

Summary of time of minima determination for V380 Cyg. 



HJD 




Type a 


Reference 


2423587.184 ±0 


.009 


P 


Kron 1935 


2423592.565 ±0 


.007 


S 


Kron 1935 


2438274.372 ±0 


.003 


P 


Semeniuk 1968 


2438279.571 ± 


.020 


s 


Semeniuk 1968 


2441256.541 ±0 


.003 


p 


Battistini et al. 1974 


2441261.620 ±0 


.004 


s 


Battistini et al. 1974 


2443791.412 ±0 


.010 


p 


Dorren & Guinan (unpub.) 


2443796.416 ±0 


.010 


s 


Dorren & Guinan (unpub.) 


2447543.936 ±0 


.005 


p 


This study (APT-PhocnixlO) 


2447548.931 ±0 


.005 


s 


This study (APT-PhocnixlO) 


2450029. 125 ±0 


.005 


p 


This study (FCC 0.8-m APT) 


2450034.033 ±0 


.005 


s 


This study (FCC 0.8-m APT) 



a P and S stand for primary and secondary minimum, respec- 
tively. 



Table 6 

Model predictions with different amounts of overshooting for the components of 

V380 Cyg. 







Primary comp. 






Secondary comp. 






Age (Myr) 


l0gT eff 


log k 2 


Age (Myr) 


logT eff 


log k 2 


0.2 


20.5 ± 1.5 


4.267 ±0.010 


-2.66 ±0.03 


20 ±3 


4.311 ±0.009 - 


-2.10 ±0.02 


0.4 


22.9 ± 1.5 


4.282 ±0.010 


-2.75 ±0.03 


23 ±3 


4.313 ±0.009 - 


-2.11 ±0.02 


0.6 


25.2 ± 1.5 


4.324 ±0.010 


-2.91 ±0.03 


26 ±3 


4.314 ±0.009 - 


-2.11 ±0.02 
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